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Abstract 

We perform a systematic X-ray analysis of six giant radio relics in four clusters of galaxies using the Suzaku 
satellite. The sample includes CIZA 2242.8-5301, Zwcl 2341.1-0000, the South-East part of Abell 3667 and previ- 
ously published results of the North-West part of Abell 3667 and Abell 3376. Especially we first observed the narrow 
(50 kpc) relic of CIZA 2242.8-5301 by Suzaku satellite, which enable us to reduce the projection effect. We report 
X-ray detections of shocks at the position of the relics in C1ZA2242. 8-5301 and Abell 3667 SE. At the position of the 
two relics in ZWCL234 1.1-0000, we do not detect shocks. From the spectroscopic temperature profiles across the 
relic, we find that the temperature profiles exhibit significant jumps across the relics for CIZA 2242.8-5301, Abell 
3376, Abell 3667NW, and Abell 3667SE. We estimated the Mach number from the X-ray temperature or pressure 
profile using the Rankine-Hugoniot jump condition and compared it with the Mach number derived from the radio 
spectral index. The resulting Mach numbers {Ai = 1.5-3) are almost consistent with each other, while the Mach 
number of CIZA2242 derived from the X-ray data tends to be lower than that of the radio observation. These results 
indicate that the giant radio relics in merging clusters are related to the shock structure, as suggested by previous 
studies of individual clusters. 
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1. Introduction 

The formation theory of galaxy clusters and groups is now 
well established: the high density peaks that exist within the 
primordial matter distribution later grow to groups and clus- 
ters by both accretion of surrounding matter and mergers of 
clumps. In the growth history of haloes, major mergers have 
a significant effect on the internal structure of clusters. X-ray 
observations have revealed a lot about the merging phenomena 
of haloes from the massive clusters scale (e.g. Markevich et al. 
2002) to galaxy groups (e.g. Kawahara et al. 2011) mainly 
by studying X-ray morphology. Shock structures induced by 
mergers were clearly found in the Bullet cluster using Chandra 
(Markevich et al. 2002). In general, however, an identification 
of a shock structure in X-ray images is difficult, except in major 
mergers including 1E0657-56, A520, A2146 and A754 (Clowe 
et al. 2006; Markevitch et al. 2005; Russell et al. 2010; Macario 
etal. 2011). 

Recent progress in radio astronomy has revealed interesting 
arc structures in merging clusters, called radio relics (Ferrari et 
al. 2008). Radio relics are believed to be a tracer of shock struc- 
tures through synchrotron radiation. Finoguenov et al. (2010) 
used XMM-Newton data of the well-known radio relic cluster 
Abell 3667 and found a sharp temperature decrement across 
the radio relic. Akamatsu et al. (2012) confirmed their result 
using data from the Suzaku satellite. 

Here, we present a systematic X-ray analysis of 6 radio 
relics in 4 clusters observed using Suzaku, including a new 
observation of a strong radio relic cluster, CIZA 2242.8- 
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5301 (hereafter CIZA 2242), which was recently discovered 
in a Giant Metrewave Radio Telescope (GMRT) observation 
(Swarup 1991). New analyses of archival data of Zwcl 2341.1- 
0000 (hereafter Zwcl2341), Abell3667SE, and a compilation 
of previous results of Abell3667NE (Akamatsu et al. 2012) and 
Abell3376 (Akamatsu et al. 201 lb) are also discussed. 

We use Ho = 70 km s"' Mpc"', Qm = 0.27 and Qa = 0.73, 
respectively. We use solar abundances as given by Anders & 
Grevesse (1989). Unless otherwise stated, the errors corre- 
spond to the 68% confidence level for a single parameter The 
rest of this paper is organized as follows. In §2, we describe the 
target selection, data reduction, and the X-ray spectral analysis. 
In §3 we derive the temperature across the relic to estimate the 
Mach number In §4 we summarize the results. 

2. Observations & Data preparation 

2.1. Targets 

In Table 3, we summarize the information about the 
relics and clusters we analyze. Four relics (CIZA2242, 
Abell3667SE, Zwcl2341N, and Zwcl2341S) are newly ana- 
lyzed in this paper and our previous analysis are used for two 
relics (Abell3667NW, A3376). We examine these clusters one 
by one. 

CIZA2242 (z = 0.192) was discovered in the Clusters in the 
Zone of Avoidance (CIZA) survey (Kocevski et al. 2007). van 
Weeren et al. (2010) reported a Mpc scale radio relic located in 
the Northern outskirts, at a distance of 1 .2 Mpc from the clus- 
ter center The radio relic is extremely narrow with a width of 
55 kpc. The spectral index at the front of relic is -0.6 + 0.05, 
which corresponds to a Mach number of Al = 4.5^^^. Although 
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Table 1. Basic properties of the clusters 





redshift 


physical length for 1' 


Temperature 


A'h 


Radio references 






(kpc) 


(keV) 


(10-" cm--) 




CIZA2242 


0.192 


190 


8.4 


334 


van Weeren et al. 2010 


Abell3376* 


0.046 


54 


4.2 


5.8 


Bagchi et al. 2006, Kale et al. 2012 


Abell3667NWt 


0.056 


66 


6.0 


4.7 


Rotttgering et al. 1997 


Abell3667SE 


0.056 


66 


6.0 


4.7 


Rotttgering et al. 1997 


Zwcl2341N 


0.270 


246 


5.7 


34 


Bagchi et al. 2002; van Weeren et al. 2009 


Zwcl2341S 


0.270 


246 


5.7 


34 


Bagchi et al. 2002; van Weeren et al. 2009 



*:Akamatsu et al. 2012 
t:Akamatsu et al. 201 lb 



Table 2. Observation log and exposure 


time after data screening 


witliout and with COR-cut(COR2 


>6GV) 


Name (ObsID) 


(R.A., DEC ) 


Observation Start time 


Exp time^ 


CIZA2242(806001010) 
CIZA2242 OFFSET( 8060020 10) 
Abell3667SE(805036010) 
Zwcl2341 (803001010) 


(340.74, +53.16) 
(339.29, +52.67) 
(303.44, -57.03) 
(355.90, +00.34) 


2011-07-28 
2011-07-30 
2010-04-12 
2008-06-27 


102.1 
52.4 
45.0 
41.0 



t:COR2 > 6GV 



this shock wave is the strongest of all the shocks in this sam- 
ple, only X-ray data from ROSAT have been available so far. 
We obtained a 120 ksec Suzaku observation of CIZA2242 per- 
formed on September 28, 201 1 with an additional 60 ksec off- 
set observation. 

A3376 (z - 0.046) is considered to be experiencing a binary 
subcluster merger. Previous BeppoSAX/PDS observations 
showed a 2.7 cr detection of a hard X-ray signal Nevalainen 
et al. 2004. However, Suzaku HXD observations gave an up- 
per limit, which did not exclude the BeppoSAX flux (Kawano 
et al. 2009). The detection of a hard X-ray signal is still con- 
troversial. Another striking feature of A3376 is a pair of Mpc- 
scale radio relics (Bagchi et al. 2006; Kale et al. 2012). Using 
Suzaku data, Akamatsu et al. (2011b) confirmed evidence for 
the presence of shock fronts across the west radio relic. 

A3667 (z - 0.0556) is well known as a merging cluster with 
irregular X-ray morphology and has large two extended ra- 
dio reUcs at the North- West (A3667NW) and the South-East 
(A3667SE) (Rotttgering et al. 1997). The overall spectral in- 
dex of North- West relic gradually varies -1.1 to -1.5 toward 
the southeast direction. A3667 exhibits an elongated X-ray 
shape to the North- West direction with an average temperature 
of 7.0 + 0.5 keV (Knopp et al. 1996; Briel et al. 2004). Recent 
XMM-Newton and Suzaku observations of the North- West ra- 
dio relic revealed a significant jump in temperature and surface 
brightness across the relic (Finoguenov et al. 2010; Akamatsu 
et al. 2012). The observed ICM properties indicate the pres- 
ence of a shock front with Mach number Ai ~ 2.4. The South- 
East radio relic has not yet been studied in X-rays. 

Zwcl2341 (z - 0.270) has two Mpc-scale structures of dif- 
fuse radio emission, which were discovered in 1 .4 GHz New 
VLA All the Sky Survey (NVSS) (Bagchi et al. 2002) at the 
North (Zwcl2341N) and South (Zwcl2341S) of the clusters. 
These radio relics are located at a Northern (Zwcl2341N ~860 
kpc ) and Southern (Zwcl2341S~ 1100 kpc) direction from 
the cluster center (Bagchi et al. 2002). Recently, GMRT ob- 



servations revealed that those radio emissions exhibit "arc"- 
like shapes (van Weeren et al. 2009). The spectral indices of 
the relics are -0.49+0.18 for Zwcl2341N and -0.76+0.17 for 
Zwcl2341S corresponding to Mach numbers Al > 3.57 and 
2.95 + 1 .39, respectively. Previous Chandra and XMM-Newton 
observations (van Weeren et al. 2009) show that the X-ray 
emission extends over ~ 3.3 Mpc in the North-South direction. 
Due to its low X-ray surface brightness, the X-ray information 
around these radio relics is still limited. 

2.2. Data reduction 

All SUZAKU observations of CIZA2242, CIZA2242 
OFFSET, Abell3667SE and Zwcl2341 were performed with 
the XIS in either normal 5x5 or 3 x 3 clocking mode. The 
detailed information about these observations is summarized 
in Table 2. For all data, 3 out of the 4 CCD chips were avail- 
able: XISO, XISl and XIS3. The XISl is the back-illuminated 
chip with high sensitivity in the soft X-ray energy range. 

We use HEAsoft version 6.11 and the calibration database 
(CALDB) dated 2011-06-30. We perform the event screen- 
ing with the cut-off rigidity (COR) > 6 GV to increase the 
signal-to-noise ratio. For CIZA2242, we set the Earth rim 
ELEVATION > 10° to avoid contamination by scattered so- 
lar X-rays from the day Earth limb. Furthermore, we ap- 
ply additional processing for XISl to reduce NXB level, 
which increased after changing the amount of charge injection 
(06/2010). The detailed processing procedure is the same as 
the one described in XIS analysis topics Regions affected 
by the calibration sources were masked out using the calmask 
CALDB file. We used data in the energy range between 0.5- 
10 keV for the FI detectors and between 0.5-8 keV for the BI 
detector For CIZA2242, we used the energy range between 
0.75-10 keV for the FI detectors and 0.6-8 keV for the BI de- 
tector. 

' http://www.astro.isas.jaxa.jp/suzaku/analysis/xis/xisI_ci_6_nxb/ 



No. ] 
(a)CIZA2242 



(b)A3376 





0.043 0.17 



.00075 0.082 0.33 0.75 



(c)A3667 



(d)Zwcl2341 





Fig. 1. X-ray image of (a)CIZA2242, (b)A3376, (c)A3667 and (d)Zwcl2341 in the energy band 0.5-8.0 keV, after subtraction of tlie NXB without 
vignetting coiTection and after smoothing by a 2-dimensional Gaussian with cr = 15 pixel =17". White and Yellow annuli indicate the spectrum analysis 
regions and the radio relic region. Magenta annuli in Panel (a) is used for the analysis of the region perpendicular to the merger axis. 



2. 3. X-ray Spectrum Analysis 

In this paper, we follow the same spectral fitting procedure as 
for the Suzaku observations of A3667 and A3376 (Akamatsu 
et al. 2012; Akamatsu et al. 2011b). The basic method of the 
spectral fit is described in subsection 4.1 of Akamatsu et al. 
(2012). 

2.3.1. Background and Foreground Emissions 

An accurate measurement of the background components is 
crucial for the ICM study in the cluster outskirts. We consider 
three sky background components: the local hot bubble (LHB), 
the Milky way halo (MWH), and the cosmic X-ray background 
(CXB). We analysed the off^set region for CIZA22420FFSET 
(ID:806002010), which is one degree away from the main 
target. For the A3667SE observation, we used the back- 
ground model reported in Akamatsu et al. 2012 Table3 and for 
Zwcl2341, we used the outer region where there is negligible 
ICM emission. 



The non X-ray background (NXB) subtracted spectra were 
fitted by the summation of the LHB: (-0.1 keV), the MWH 
(~0.3 keV) and the CXB as apec + wabs(apec + powerlaw). 
We fix the redshift and abundance of all apec components to 
zero and unity, respectively. The NXB component was esti- 
mated from the dark Earth database by the xisnxbgen FTOOLS 
(Tawa et al. 2008) and was subtracted from the data before 
the spectral fit. We used the XIS response assuming uniform 
brightness on the sky by xissimarfgen (Ishisaki et al. 2007). 
The best fit parameters and spectrum of CIZA22420FFSET 
are shown in Table 3 and Fig. 3, respectively. 
2.3.2. Spectrum fitting 

We extracted spectra from regions defined by their 
distances from the center (22H2™41.'^9, 50°03™00^) 
2.'5, 5.'0, 7.'7, 10.'3, 12.'8, 16.'7 (Fig la ) for CIZA2242, 
16.'8, 21.'0, 25.'2 centered on (20M2'"3P, -56°49"12^) 
(Fig lb ) for A3667SE, 1.'5, 4.'0, 8.'5, lO.'O centered on 
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(d)7.'7-10.'3 



Energy (keV) 





(e)10.'3-12.'8 



Energy (keV) 





(f)12.'8-16.'7 



Energy (keV) 




Energy (keV) 



Energy (keV) 



Energy (keV) 



Fig. 2. The NXB subtracted spectra of six atmuli of the CIZA2242 (corresponding to white annuli in Figure 1 . The XIS BI (BlacJc) and FI (Red) spectra 
are fitted with the ICM model (wabs + apec), along with the sum of the CXB and the Galactic emission (apec + wabsfapec + powerlaw)). The ICM 
components are shown in magenta. The CXB component is shown with a black curve, and the LHB and MWH emissions are indicated by green and blue 
curves, respectively. The total background components are shown by the orange curve. 
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Fig. 3. The spectrum of the CIZA22420FFSET used for the back- 
ground estimation, after subtraction of NXB and the source. The XIS 
BI (Black) and FI (Red) spectra are fitted with CXB + Galactic compo- 
nents (LHB, MWH) (apec+wabs(apec+powerlaw)). The CXB spec- 
trum is shown with a black curve, and the LHB and MWH components 
are indicated by green and blue curves, respectively. 

(23H3'Ml\-i-00°18"'190forZwcl2341N, 1.'5, 4.'0, 7.'0 cen- 
tered on (23H3'"4P,H-00°18"19^) (Fig Id ) for Zwcl2341S, 
respectively. For CIZA2242, the center was determined by 
visually fitting a circle to the radio relics. For other clusters, 
the center was determined from their X-ray surface brightness 
peaks. 

The response matrix files (RMFs) and ancillary response 
files (ARFs) were generated using xisrmfgen and xissimarf- 



Table 3. Summary 
ters for each 


of the 
background 


best fit 
observations. 


parame- 




C1ZA2242 


A3667SE 


Zwcl2341 


LHB 

ytr(keV) 
norm* (xlO"'') 


0.08 (fix) 
12 26 " 


O.lliO.Ol 
6.85!!™ 


0.10+0.01 

4 85+' ™ 


MWH 
fcr(keV) 
norm* (xlO"'*) 


0.29+0.01 
55.70!^-7 


0.37 + 0.03 


0.25+0.02 

1 l+"-3 


CXB (r= 1.41) 

«orm"'" (xlO"'*) 


8.8 +0.4 


9.6+0.6 


11.2+0.3 



*:Normalization of the apec component scaled with a factor l/400;r. 
Norm=TO / ncnHdV/(4;r(l +z2)D2 )x IQ-''' cm-^ arcmin-^, 
where is the angular diameter distance to the source. 
The refen'ed CXB intensity normalization in Kushino et al. 2002 is 
9.6 X lO"-* for F = 1.41 in the unit of photons keV"' cm"^ s"' at 1 keV. 



gen (Ishisaki et al. 2007). Since the Galactic and CXB com- 
ponents can be regarded as being spatially uniform in the XIS 
field of view, we use an uniform ARF for these components. 
For the ICM emission, we need to put the flux distribution into 
xissimarfgen. For A3667 and A3376, we adopted the j6-model 
distribution with y6=0.54, rc=2.'97 and p = 0.40, = 2.'03, re- 
spectively. For C1ZA2242 and Zwcl2341, we used the Suzaku 
XIS image (0.5-8.0 keV). 

The spectral fitting of each region was performed as de- 
scribed below. We modeled the NXB subtracted spectrum as 
the sum of the ICM emission and each background component 
discussed in Sec 2.3.1. Intensity, temperature, photon indices 
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of CXB and galactic components were all fixed at the value de- 
termined in offset observation. The interstellar absorption was 
kept fixed using the twenty-one centimeter line measurement 
of the hydrogen column, A^h = 33.4, 4.7 and 3.4 x 10^° cm"^ 
for CIZA2242, A3667SE and Zwcl2341, respectively (Dickey 
& Lockman 1990). We simultaneously fitted both the BI and 
FI spectra with XSPEC verl2.7. For the inner regions of the 
radio relic, the free parameters were temperature kT, normal- 
ization norm and metal abundance Z of the ICM component. 
Otherwise, we fixed the ICM metal abundance to the typical 
value in cluster outskirts, Z = O.2Z0 (e.g. Fujita et al. 2008). 
The relative normalization between the two sensors was a free 
parameter in this fit to compensate for cross-calibration errors. 

3. Results 

3.1. ICM properties of individual cluster 

3.1.1. CIZA2242 

We first perform the spectral fitting within the radio relic (~ 
1.2 Mpc= 6'5): the temperature and the metal abundance are 
kT = 7.88 + 0.20 keV and Z = 0.28 ± 0.03, respectively. The 
unabsorbed flux and luminosity within the radio relic in the 
0.5-10 keV band are Fo.5-10.0 keV - 8.5 x 10"'~ erg s cm"- and 
i'O.s-io.o keV = 8.4 X lO'*'* erg s, respectively. 

Then, we performed the spectral fitting for each annulus ex- 
amined in Sec 2.3.2. Figure 2 displays spectra of six annuli. 
Figure 4(a) shows the radial profile of the ICM temperature 
in CIZA2242. The temperature within 6. '5 (1.25 Mpc) from 
the cluster center shows a fairly constant value of ~ 8 - 9 keV. 
As mentioned in section 1, van Weeren et al. (2010) reported 
the sharp radio relic at that radius, which is indicated by ver- 
tical dashed line in Panel (a). We found that the temperature 
profile exhibits a significant jump from 8.3 keV to 2.1 keV 
across the radio relic. Because the ICM emission is lower than 
the background level in the outer bins (r > 7.7'), we consider 
the variability of the CXB as a systematic error The fluctu- 
ation of CXB intensity was estimated based on our previous 
work (Akamatsu et al. 201 la). The estimated fluctuations span 
15-25 %. The resultant parameters after taking into account the 
systematic error are shown blue dotted lines in figure 4 (a). The 
temperature jump is consistent with the prediction of numeri- 
cal simulation by van Weeren et al. (201 1). For a reference, we 
plotted the "universal" temperature profile expected from the 
scaled temperature profiles (Burns et al. 2010). 

As shown in Fig. 5(a), we could not see the surface bright- 
ness jump across the relic directly in the profile. However, in 
Appendix we show that this is caused by the limited Suzaku 
spatial resolution, and the observed profile is in fact consistent 
with the expected density jump 

3.1.2. A3376 

Previous studies of A3376 found that the global mean tem- 
perature is 4.0 keV and showed the presence of a pair of Mpc- 
scale radio relics consistent with 1 .4 GHz VLA NVSS observa- 
tions (Bagchi et al. 2006). A detailed study of the A3376 radio 
relic has already been published by Akamatsu et al. (2011b). 
We refer to that paper for the images and X-ray information 
and we will include A3376 in the discussion. 



3.1.3. Abell3667NW 

Suzaku performed 3 pointing observations of Abell 3667 
along the North- West merger axis in May 2006. Figure. 1(c) 
displays the NXB subtracted X-ray image of A3667. A detailed 
study of A3667NW radio relic has already been published by 
Akamatsu et al. (2012). We refer to that paper for the images 
and X-ray information and we will include A3667NW in the 
discussion. 

3.1.4. A3667SE 

Figure. 4(d) indicates the radial profile of the ICM tempera- 
ture in the South-East direction of A3667, respectively. The in- 
nermost 2 bins of the cluster are identical to that of A3667NW. 
The temperature within the SE radio reUc (~ 1 .3 Mpc) from the 
cluster center shows a fairly constant value of ~ 7 keV, which is 
consistent with the NW direction (Akamatsu et al. 2012). Our 
temperature profile exhibits a jump across the SE radio relic 
region. Although we could not see the surface brightness jump 
(Fig. 5d), the pressure profiles (Fig. 6d) exhibit a slight drop 
due to the temperature jump. 

5.7.5. Zwcl2341N&Zwcl2341S 

We first extract the spectrum within 3' (~ 740 kpc) and fix 
the redshift value to 0.270. The spectrum was well fitted with 
the above model (xl - 1.04 for 473 degrees of freedom). We 
find a temperature and metal abundance of kT = 5.16 + 0.32 
keV and Z = 0.25 + 0.08, respectively. The unabsorbed flux 
and luminosity in the 0.5 - 10 keV band are /^o.s-io.o keV = 1 .0 x 
10"'^ erg s cm"^ and Lo.5-10.0 keV = 2.2x lO'*'* erg s, respectively. 
The derived temperature is consistent with previous Chandra 
observations (~ 5 keV) (van Weeren et al. 2009). 

The radial profiles of the ICM temperature are shown in 
Fig. 4(e) and (f) for Zwcl2341N and Zwcl2341S, respec- 
tively. The temperature drops from 4keV to 1 keV toward 
the Northern direction or from 4 keV to 2 keV toward the 
Southern direction. Figure 5(e,f) and Figure 6(e,f) show the 
surface brightness and the pressure profiles of Zwcl2341N and 
Zwcl2341S. We could not see any jump in our wide radial bins 
for temperature, surface brightness, and pressure profiles for 
this cluster 

3.2. The scaled temperature profiles 

We normalized the radial temperature profile of each clus- 
ter by the flux- weighted average temperature. The resulting 
scaled temperature profiles are shown in Fig 7. We find that 
the temperature profiles of four relics, including the previous 
two results, exhibit significant discontinuities across the ra- 
dio relic. These discontinuities can be associated with the 
shock front found in the radio relics. For the other two sam- 
ples (Zwcl2341N,S), we could not see any jump in temperature 
profiles. The angular resolution of Suzaku was not enough to 
check the presence of a discontinuity. Therefore, we focus on 
these four examples of temperature discontinuity. 

Another remarkable point of the scaled temperature pro- 
file is its flatness within the radio relic. Recent X-ray results 
show that the ICM in relaxed clusters has the "universal" de- 
cline temperature profile up to virial radius (Markevitch et al. 
1998; De Grandi & Molendi 2002; Vikhlinin et al. 2005; Pratt 
et al. 2007; Akamatsu et al. 2011a). The flatness of temper- 
ature within the radio relic might need a heating mechanism, 
such as shock heating. 
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Fig. 4. X-ray temperature profiles for the six radio relics. Vertical dashed lines indicate the position the of radio relics. Gray dotted lines indicate the 
"universal" temperature profile expected from the scaled temperature profiles (Bums et al. 2010). The red and blue horizontal bars show the pre and post 
shock quantities we use to derive Mach number (Sec 3.3). The vertical bars indicate a Icr error. In (a)CIZA2242, the gray crosses show the data for a 
reference sector as shown in Fig. 1(a) magenta annuli. The uncertainty range due to the fluctuations of CXB intensity are shown by two blue dotted lines. 



Table 4. The post- and pre- shock ICM quantities and Mach number derived from X-ray and radio observations. 



paste shock 

(keV) 



pre shock 

(keV) 



Mx, 



kT 



>^radio 



3.15+0.52!'i'-^!!* 
2.94+0.60" 
2.41+0.39 
1.75+0.13 



CIZA2242 
Abell3376 
Abell3667NW 
Abell3667SE 



8.33+0.80 + 0.40* 
4.81+0.29 
5.52+0.95 
6.34+0.38 



1.35. 
2.03. 
3.59= 



:0.35 
:0.34 
:0.28 



-0.60+0.05 4.58+1.32 

-1.00+0.10 2.23+0.20 

-1.10+0.20 2.08+0.37 

-1.50+0.17 1.73+0.58 



*: The first and the second errors are statistical and systematic (CXB fluctuations). 



As mentioned above, the temperature of relaxed clusters 
tends to have the "universal" profile. It is natural that the flat 
temperature profiles of merging cluster will settle to the uni- 
versal profile in less than a Hubble time. In our samples, the 
locations of the radio relic are far from cluster center, typically 
rshock ~ Mpc (0.5r2oo). Due to the low density («(. = lO"*"^ cm"-') 
in the outskirts of clusters, the cooling time in those regions ex- 
ceeds the Hubble time. We might consider other processes to 
achieve the universal profile. 

3.3. Mach number at Discontinuity 

We estimate the Mach numbers using temperature jumps 
across the relics. The Mach number can be obtained by apply- 
ing the Rankine-Hugoniot jump condition, assuming the ratio 



of specific heats as y = 5/3, as 

(2) 

T\ , T2 are the pre-shock and post-shock temperatures and 
Mx,kT is the estimated Mach number. The post and pre quan- 
tities we use are marked by red and blue horizontal bars in 
Figure 4. 

We compare our results with the Mach number derived 
from previous radio observations. Based on Diffiisive Shock 
Acceleration (DSA) theory, we estimate the Mach number 
based on the radio spectrum, as 
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Fig. 5. Surface brightness profiles (0.5-10 keV). Vertical dashed lines indicate the position of radio relics. 



a = - 



3M--,. +1 

radio 



(3) 



where a is radio spectral index and Alradio is the Mach number 
from the radio observation. For strong shocks, the radio spec- 
tral index will be flat about -0.5. We used radio spectral indices 
from the radio observations summarized in Table 1 . The result- 
ing Mach numbers are shown in Table 4. Figure 8 displays the 
Mach number from the X-ray temperature and radio observa- 
tions. As shown in this figure, we find that the two different 
methods give consistent results within the current observational 
noise. 

If we assume in the case of CIZA2242 that the mach num- 
ber derived from the radio observation is correct, the expected 
post shock temperature would reach kr = 15 keV. Although 
our results are well reproduced by a 1 temperature model, we 
also fit 2 temperature models to the spectrum extracted from 
the 5.'0- 7. '7 annulus, which is located in the region where 
the shock just passed. The metal abundances of the two apec 
components were tied to be the same, and only the temperature 
and normalization were set free. We then obtained a statisti- 
cally acceptable fit with^^/ d.o.f. - 228/244 compared with the 
single temperature case of 231/246. The rather low statics of 
the spectrum hampered us to distinguish between the 2kT and 



IkT models. The hot component temperature was derived as 
fcrhigh = 1 1 .5+^-^^ keV, and the cool one as /tr,ow =2.05+^-^] keV. 
The intensity of the cool component is about a order lower than 
hot one. Based on those analyses, there were no signatures of 
the presence of high-temperature components. 

4. Summary 

From simulations of the structure formation of the 
Universe (Enfilin & Briiggen 2002; Vazza et al. 2009), large 
accretion shocks and low merger shocks around clusters of 
galaxies are expected. To understand the physics of cluster 
evolution, cluster merger events are important process. Suzaku 
XIS is best suited instrument for measuring the ICM emission 
beyond the radio relics thanks to its low and stable detector 
background. In this paper, we analyse six radio relics of four 
clusters with Suzaku XIS and measured their temperature pro- 
files. We found clear temperature jumps across four relics and 
derived the Mach number from the pre- and post temperature. 
The derived Mach number is almost consistent with that from 
the radio observation. These results strongly suggest that a 
shock structure exists at the radio relic. 

We are grateful to Takaya Ohashi and Jelle de Plaaa for in- 
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Fig. 6. Pressure profiles derived by the deprojection technique as described in Akamatsu et al. (2011a). Veilical dashed lines indicate the position of 
radio relics. The vertical bars are 1 cr errors. 
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Japan Society for the Promotion of Science (JSPS) Fellows 
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Appendix 1. Surface brightness jump in CIZA2242 

Because of the non-detection of surface brightness jump in 
CIZA2242, we performed an additional check. Based on previ- 
ous works (e.g., Markevich et al. 2002), we modeled the elec- 
tron density profile switching at radio relic (r = 6.'5) from the 
beta model to a power law. In this model, we assume that two 
cases of the jump amplitude, corresponding to the Mach num- 
ber Al = 2.0 (our lower limit including the systematic uncer- 
tainty, see Table 4) and M=4.5 (the radio value). Assuming 
spherical symmetry, we calculated the square of density as sur- 
face brightness and convolved it with the Suzaku PSF. The 
profile was binned with the same size as the observed sur- 
face brightness. Because the jump in the surface brightness 
is important, we selected the profile across the radio relic 
(r - 5.'0 - 10.'3). Finally, we normalized the profile to the 
observed value at the post shock region (r = 5.'0 - 7.'7). The 
resulting profiles are shown in Fig. 9b (Blue and Black). The 
observed brightness profile, shown by red crosses, is consistent 
with the interval of M obtained from our temperature measure- 
ments. 

For the CIZA2242 case, the PSF of Suzaku corresponds to 
380 kpc, which is probably much larger than the length of 
shock. Hence the surface brightness jump should be signifi- 
cantly diluted by other area in the bin. We need more obser- 
vations with a higher angular resolution to confirm the shock 
structure in the surface brightness profile. 
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